Shot noise in semiclassical chaotic cavities 
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We construct a trajectory-based semiclassical theory of shot noise in clean chaotic cavities. In 
the universal regime of vanishing Ehrenfest time te, we reproduce the random matrix theory result, 
and show that the Fano factor is exponentially suppressed as te increases. We demonstrate how our 
theory preserves the unitarity of the scattering matrix even in the regime of finite te. We discuss 
the range of validity of our semiclassical approach and point out subtleties relevant to the recent 
semiclassical treatment of shot noise in the universal regime by Braun et al. cond-mat/0511292 . 
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Introduction. Quantum transport through chaotic 
ballistic cavities is often well described by Random Ma- 
trix Theory (RMT) 0. Despite its many successes, or 
should we say, because of these successes, one might won- 
der what is the origin of this RMT universality, and under 
what conditions do system specificities modify the RMT 
of transport. System specific contributions to transport 
originate from the underlying classical dynamics, which 
suggests that one employs semiclassical methods based 
on classical trajectories Q. Indeed, the semiclassical pro- 
gram toward a microscopic foundation for the RMT of 
transport, including explicit bounds for its regime of ap- 
plicability, is currently on its way to being completed 
successfully @ S H S 0, S @ 

Here we contribute to this program by deriving the 
zero-frequency shot noise power S for quantum chaotic 
systems. The interest in shot-noise, the intrinsically 
quantum part of the fluctuations of a non-equilibrium 
electronic current, is that it often contains information 
on the system that cannot be obtained through con- 
ductance measurements. For instance, shot-noise exper- 
iments have determined the charge and statistics of the 
charge carriers in superconducting heterostructures or in 
the fractional quantum hall effect 01 • I n this paper, 
we consider an open ballistic quantum dot |llj | carry- 
ing a large number of conducting channels and accord- 
ingly neglect electron-electron interactions. We repro- 
duce the RMT result, and show how shot noise deviates 
from RMT predictions in the semiclassical limit. We 
calculate the Fano factor F — S/S p , given by the ra- 
tio of S to the Poissonian noise S p = 2e(J) that would 
be generated by a current flow of uncorrelated electrons. 
According to the scattering theory of transport one has 
F = Tr[ttt(l - ttt)]/Tr[t+t] 03. If one makes the RMT 
assumption that the transmission matrix t is the JVl x 
off-diagonal block of a (N^+N-n) x (Nl+N-r) random uni- 
tary scattering matrix, one gets F = Nl,Nr/(Ni, + Nr.) 2 
0, liol ] , in term of the number of quantum channels JVl 
and Nn carried by the contacts to the left and right 
leads. Ref. p| carried out the first semiclassical calcu- 
lation of F for the specific case of quantum graphs. The 
difficulty is to calculate Tr[ttttH] = T,i,j, q l*j\i| 2 IV*| 2 + 
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FIG. 1: (Color online) The four dominant contributions to 
Tr^tt^t]. Paths are paired everywhere except at encounters 
where two of them (71, 73) cross each other (solid lines) while 
the other two (72,74) avoid the crossing (dashed lines), (a) 
Contribution D\ has uncorrelated escape on both sides of the 
encounter, (b) Contribution D2 and D3 have correlated es- 
cape only on one side of the encounter, (c) Contribution D4 
has correlated escape on both sides. 



Eij.J^il^i.pl' + E^i^^^p^i- Ref. em- 
ployed a diagonal approximation to calculate the first 
two terms and identified the dominant four-trajectory 
contributions to the third one. Quantum graphs funda- 
mentally differ from continuum models which we treat 
here. In our semiclassical derivation we find the domi- 
nant contributions to F from the path pairings shown in 
Fig. ^ These pairings are similar to those considered in 
Ref. [f| for quantum graphs. However, unlike quantum 
graphs, chaotic systems have continuous families of scat- 
tering trajectories with similar actions, which means in 
particular that we cannot make a diagonal approximation 
to evaluate the contributions F>i and D3 shown in Fig. ^ 
This important point was not addressed in Ref. . 

Exploring the range of validity of RMT for chaotic sys- 
tems, we find F to be exponentially reduced [T^| . 



F = N L N K (N h + N R )~ 2 exp[-T£/T D ], 



(1) 



2 



for systems with left (right) lead width, Wl (Wr), such 
that the width of leads Wl,r ^l^L. These systems 
witness the emergence of the new Ehrenfest time scale 
rg = A -1 lnf?!^ (t{/td) 2 ], which generically induces sig- 
nificant deviations from the RMT of transport . Here, 
^cff = %/ (pfL), L is the linear system size, pp the Fermi 
momentum of the particle with mass m, Tf the time of 
flight, rrj the dwell time through the system, and A the 
Lyapunov exponent of the chaotic classical dynamics. 

Our semiclassical calculation correctly captures both 
the universal regime with TjS/td <C 1 and the deep semi- 
classical regime where Tg becomes comparable to or ex- 
ceeds td . We reproduce Eq. (JTJ and explicitly show that 
the exponential suppression of F is due to paths shorter 
than which become noiseless 0, 0, EI ■ We demon- 
strate the unitarity of the theory by calculating both 
F = Tr[ttt(l - ttt)]/Tr[ttt] and F = Tr[tttrt r ]/Tr[ttt]. 
We finally comment on the current limitations of the 
trajectory-based semiclassical approach. 

We consider a two-dimensional chaotic quantum dot 
ideally connected to two external leads. We require that 
the size of the openings to the leads is much smaller 
than the perimeter of the system but is still semiclas- 
sically large, 1 <C N^,N R <c L/Xf- This ensures that 
the chaotic dynamics inside the dot has enough time to 
develop. The system's transport properties are given by 
its scattering matrix S, with an TVl X JVr transmission 
block t, and an N-^ x Nz, reflection block r. To calculate 
the Fano factor, one needs to calculate the conductance 
g = Tr[tH], as well as Tr[tHttt]. 

Semiclassically, the transmission matrix reads |l5j . 

tji = -(2^)- 1 / 2 /d 2/0 /dyV(d^/dy ) 7 /2 

x(i|2/)(yo|i)exp[iS' 7 /a + i7r/i 7 /2] , (2) 

where \i) is the transverse wavefunction of the ith lead 
mode. This expression sums over all paths 7 (with clas- 
sical action S 7 and Maslov index /i 7 ) starting at yo on 
a cross-section of the injection (L) lead and ending at y 
on the exit (R) lead. We approximate ^2 n (y'\n){n\y) w 
<>(y' — y) 0i to write Tr^tt^t] as a sum over four paths, 
7I from ?/oi to yi, 72 from y 03 to y lt j3 from y Q3 to y 3 
and 74 from yoi to y 3 , 



Tr^tt^t] 



1 



dyoidj/03 / dyidy. 



(2^) 2 JL JR 
x V] A 1 ^A iZ A 1 - 1 A 1 y exp[i£S/7i] . (3) 

7 l,"- 7 4 



1/2 

Here, A 7 = [dp^/dyo^ and SS = 5 7 i — 5 7 2 + S l3 — 5 7 4 
(we absorbed all Maslov indices into the actions 5 7 ,). 
We are interested in quantities averaged over variations 
in the energy or the system shape. For most contribu- 
tions, SS/h oscillates wildly with these variations. The 
dominant contributions that survive averaging are those 
for which the fluctuations of 5S/h are minimal. They 
are shown in Fig. ^ Their paths arc in pairs almost 



everywhere except in the vicinity of encounters. Going 
through an encounter, two of the four paths cross each 
other, while the other two avoid the crossing. They re- 
main in pairs, though the pairing switches, e.g. from 
(7l;74) and (72573) to (71572) and (73574) in Fig.QJ,. 
Paths are always close enough to their partner that their 
stability is the same. Thus, for all pairings in Fig. ^ 
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73^7! ■ 



(4) 



We define P(Y, Yo; t)6yS95t as the product of the mo- 
mentum along the injection lead, pf cos#o, and the clas- 
sical probability to go from an initial position and angle 
Yo = (yo, Oq) to within (Sy, 69) of Y in a time within St 
of t. Then the sum over all paths 7 from y to y is 



'dtjd9 jd9 P(Y,Y ;i) [• 



r -(5) 



For an individual system, P has (5-functions for all classi- 
cal trajectories. However averaging over an ensemble of 
systems or over energy gives a smooth function 



(P(Y,Y ;f)) 



Pf cos 9 cos 9 
2{W L + W r )td 



exp[-t/Ta]. (6) 



Using Eqs. JHJ) and (jSJ to calculate the conductance 
within the diagonal approximation directly leads to the 
Drude conductance (Tr[ttt]) ~ g D = N L N R /(N L + N R ). 
This level of approximation for (Tr[ttt]) is sufficient to 
obtain F to leading order in N^^. We now use Eqs. J3J), 
(0} and JSJ to analyze the contributions in Fig. 2] 

There are two things that can happen to two pairs of 
paths as they leave an encounter. The first is uncorrelated 
escape. The pairs of paths escape when the perpendic- 
ular distance between them is larger than Wl,r, which 
requires a minimal time 7V( e )/2 = A -1 ln[e _1 iy/L] be- 
tween encounter and escape. The two pairs of paths 
then escape in an uncorrelated manner, typically at com- 
pletely different times, with completely different mo- 
menta (and possibly through different leads). The second 
is correlated escape. Pairs of paths escape when the dis- 
tance between them is less than Wl,r, then the two pairs 
of paths escape together, at the same time through the 
same lead. 

Contributions to the Fano factor. Taking into 
account the two escape scenarios just described, we write 
(Tr[tttttt]) = D 1 + D 2 + D 3 + D A . Each of these four 
contributions, sketched in Fig. ^ can be written as 



(2irh) 2 



dY i dY 03 /dYi dY 3 /dt x dt 3 

L JR 



x (P(Y 1 ,Y i;ti) P(Y 3 ,Y 03 ;t3)> exp[iSS Di /h] , (7) 

where subscripts 1,3 make the connection to Fig. ^ 
When evaluating Eq. Q the joint exit probability for 
two crossing paths has to be computed. 

To evaluate D\, we use the method developed by 
Richter and Sieber Q, while taking into account that 
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paths in the same region of phase-space (shaded areas in 
Fig. have highly correlated escape probabilities • 
Here the action difference is (JSz^ = Epe 2 /X Q, where e 
is the crossing angle shown in Fig. [TJi. We write 

■P(Yj, Yo»; t») = / dRi-P(Yi, R,; ij — i£)P(Rj, Yoi;4) , 

where P is the probability for the classical path to exist 
(not multiplied by the injection momentum), and R; is a 
point in the system's phase-space (r^,^) visited at time 
t\, with <pi giving the direction of the momentum. Wc 
choose Ri and R3 as the points at which the paths cross, 
so R3 = (n, <pi ± e) and dR3 = v F sin edt[dt 3 de. Thus 

Di = 2(2^)- 2 ^dY idY 03 

x rdeRc[c i5s - 1 /' i ](/(Yo 1 ,Yo3;e)). (8) 
Jo 

/(Y01, Y03; e) is related to the probability that 73 crosses 
7I at angle ±e. Its average is independent of Yoi.,03, so 
(J(Y i, Y03; e)) = (1(e)). For Di, injections/escapes are 
more than 7V(e)/2 from the crossing, so 



(/(e)) = 2vl sine J dYidY 3 J dRi 



rti-T/2 rtx, rt 3 -T/2 

/ dti / dti / dt 3 / di 3 

/T JT/2 JT JT/2 

x (P(Yi,Ri;t 1 -ti)P(Ri,Y 01 ;ti) 

x P(Y 3 ,R 3 ;t 3 -t' 3 )P(R 3 ,Y 03 ;t' 3 )), (9) 

where T is shorthand for IV (e). We next note that 
within 7V(e)/2 of the crossing, paths 7I and 73 are so 
close to each other that their joint escape probability is 
the same as for a single path (this was absent from Ref. Q 
and was first noted in Ref. [17| )• Elsewhere 71, 73 escape 
independently through cither lead at anytime, hence 



Pf t d N r cos #01 cos £>03 sin e 



-7V(e)/TD 



irhm 



{N L + 7V R )3 



(10) 



where we used N R — (ttK) VfWr, and assumed that 
the probability that 73 is at R3 at time t' 3 in a system 



of area A is (2nA)~ 



i[2irhT D (N L + Nr)}- 1 . Then 



the Yoi^s-integral in Eq. © gives (2IF L ) 2 , while the e- 
integral is dominated by e < 1 and yields a factor of 
-7r^(2 J B F r D )- 1 e- T EAD| 1 + C ,[( ArD )-i]} g_ Thug 



-N 2 N 2 (N L + N R )~ 3 exp[-r|/r D ] . (11) 



The contribution D2 is shown in Fig. with Fig. [21 
showing the paths in the correlated region in more detail. 
Noting that 72 decays exponentially towards 7I, we find 
the action difference between the two paths to be 

S2-S1 = pf{vqi - 2/03) sin 601 

+ \mX(y 01 - ymf cos 2 6 0X . (12) 




FIG. 2: (Color online) Paths for the D2 and D4 contribu- 
tions when they are in the correlated region (close to L lead). 
Paths 7I and 73 (solid black lines) start on the cross-section 
of the L lead at positions 2/01 and yo3 with transverse mo- 
menta pp sin #01 and pf sin #03, respectively. In the basis par- 
allel/perpendicular to 7I, the initial position and momentum 
of path 73 are r ± = (2/01-2/03) cos #01, r || = (2/01 -2/03) sin 6 01 
and po± ~ Pf(#oi — #03)- Contribution D3 has exactly the 
same structure close to the R lead. 



The equation for S4 — S 3 has the opposite sign for (7/01 — 
2/03) and 6»oi replaced by 6 03 . In terms of (ro±,Pox) 

6S D . 2 = -(po± + mXr ±)r ± , (13) 

where we have dropped cubic terms (they only give h- 
corrections to the stationary-phase integral). We next 
perform the average in Eq. (JJJ. We define T{ v (r ±,p ±) 
as the time for which 7I and 73 are less than W apart, 
and insist that the paths are more than W apart before 
they escape to the right. Hence we must evaluate 



dYxdY 3 / dt 1 dt 3 (P(Y 1 ,Yoi;ti)P(Y3,Yo3;t 3 )) 



N&p F cos 6»oi cos 6*03 r „, , w , nA , 

= Tat ,at \2 exp[-T w (r 0± ,pQ ± )/T B } . (14) 

(N L + N R ) Z 

Inserting this into Eq. (|7J), we change integration vari- 
ables using pf cos 6*03 dY 03 = dr j^dp _L [Oli and then 
define po = po± + mXrox- In the regime of interest 
T{y(rox,pox) — A -1 hi[(mXW)~ 1 po}- Evaluating the in- 
tegral over ?*o^ leaves a p~o-intcgral which wc cast as Eulcr 
T-functions. To lowest order in (Atd) -1 we find, 

D 2 = iV L iV 2 (7V L + iV R )- 2 cxp[-rg/T D ]. (15) 

Substituting TVl <-> N R in the derivation of Eq. I|15|) gives, 

D 3 = NlN R (N L + N R y 2 eM-TE/Tn]- (16) 

The contribution D4 is shown in Fig. with Fig. [2] 
showing the paths in detail at the L lead. This con- 
tribution can be evaluated in a way similar to D2, 
the difference being that the paths escape before time 
TyyfroXiPax), he. before becoming a distance W apart. 
The paths are always correlated, so the escape probabil- 
ity for the two paths equals that for one. Moreover, both 
paths will automatically escape through the same lead, 
hence 

dY!dY 3 f "dt 1 dt 3 (P{Y 1 ,Y i;t 1 )P{Y 3 ,Y 03 ;t 3 )) 



N R pl cos 6>oi cos 6>Q3 
N h + N K ' 



1 



-Tl v (r ±,p ±')/Tn 



)■ (17) 
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Performing the same analysis as for D 2 we find that 

D 4 = iV L Ar R (Ar L + Ar R )-i(i_exp[-T£/T D ]). (18) 

The Fano factor is given by F = 1 -g^ 1 (D 1 +D 2 + D 3 + 
D 4 ). Our results of Eqs. {I5Jl, (EJl and {TBJ) show that 
D2 + -D3 + -D4 = gu- One hence gets F = — -Di/<?d- From 
Eq. (|11|) . one finally obtains our main result, Eq. fij. 
The splitting of phase-space discussed in Refs. [y, 0, L1J] 
for Tg^To naturally emerges here. For paths shorter than 
Tg, only D4 is non-zero. This cancels these path's Tr[t^t]- 
contribution, making them noiseless. 

Preservation of Unitarity. The unitarity of the 
scattering matrix ensures that t't + r^r = 1 and hence 
the Fano factor can be written as F = g^ 1 (Tr[t+trt r]) . 
We calculate this expression to explicitly show that our 
method preserves unitarity. We first note that there is 
no contribution D3 nor D 4 to Tr[t+trtr]. We are left 
with the calculation of two contributions, D[ and D' 2 , ob- 
tained from Di and D 2 shown in Fig. ^i,b with y 01 , y 03 
and 2/3 on the left lead and y\ on the right lead. The 
calculation proceeds as for D\ and D 2 , with one factor 
of N R /(N L + N R ) replaced by N L /(N L + N R ) in both 
contributions. The sum of these two contributions is 
D[ +D' 2 = e- T v/ T °NlN R {N L + N R y 3 , the Fano factor 
is then F = (D[ + D' 2 )/gu, which reproduces Eq. (J). 

Off-diagonal nature of all contributions. In our 
analysis we allow for the fact that open chaotic systems 
have continuous families of paths with highly correlated 
actions coupling to multiple lead modes. For example 
paths 7I and 72 in Fig. [21 have an action difference 
given in Eq. (|12fl . which does not fluctuate under en- 
ergy or sample averaging. The stationary phase integral 



for £>2,3,4 over such paths is dominated by paths 7I and 
73 with po± ~ — m\ro±. Since ro± is integrated over the 
width W of the lead, po± varies over a range of order 
mXW, these contributions are clearly not diagonal in the 
lead mode basis. Upon completion of this manuscript, we 
became aware of Ref. Q, which presents a semiclassical 
calculation of F for = 0. While their method is super- 
ficially similar to ours, they make a diagonal assumption 
to get the contributions that we call £>2,3,4- Our analy- 
sis shows that this is unjustifiable. Such an assumption 
would moreover violate unitarity for finite rg. 

Regime of applicability of these semiclassics. 
We appear to be the first to report that all trajectory- 
based semiclassical methods used so far in the theory of 

transport (including in the present article) arc only appli- 

1/2 

cable in the regime W > h cS L. Dominant off-diagonal 

contributions such as those discussed above have encoun- 

1/2 1/2 
ters of a typical size ~ fij^ L. When W < TiJ^ L, the two 

non-crossing paths (i.e. 72 and 74 in Fig. at an en- 
counter are a distance apart greater than W . The prob- 
ability that one of the four paths escapes while the other 
three paths remain in the system is of order te/td, where 
te ~ A" 1 ln[7i ff (L/W) 2 ] is the time over which this path 
is a distance of order W from any of the other paths. The 
current methods fail once this is taken into account, sug- 
gesting that diffraction effects may become important. 
We believe that the regime 7j e ff < (W/L) < fij^ is well 
described by RMT, and thus suspect this diffraction may 
be the microscopic source of RMT universality in this 
regime. Clearly this regime merits further study. 
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